Pulsars Cannot Account for the Inner Galaxy's GeV Excess 
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Using data from the Fermi Gamma-Ray Space Telescope, a spatially extended component of 
gamma rays has been identified from the direction of the Galactic Center, peaking at energies 
of ~2-3 GeV. More recently, it has been shown that this signal is not confined to the innermost 
hundreds of parsecs of the Galaxy, but instead extends to at least ~3 kpc from the Galactic Center. 
While the spectrum, intensity, and angular distribution of this signal is in good agreement with 
predictions from annihilating dark matter, it has also been suggested that a population of unresolved 
millisecond pulsars could be responsible for this excess GeV emission from the Inner Galaxy. In this 
paper, we consider this later possibility in detail. Comparing the observed spectral shape of the Inner 
Galaxy's GeV excess to the spectrum measured from 37 millisecond pulsars by Fermi, we find that 
these sources exhibit a spectral shape that is much too soft at sub-GeV energies to accommodate 
this signal. We also construct population models to describe the spatial distribution and luminosity 
function of the Milky Way's millisecond pulsars. After taking into account constraints from the 
observed distribution of Fermi sources (including both sources known to be millisecond pulsars, 
and unidentified sources which could be pulsars), we find that millisecond pulsars can account for 
no more than ~10% of the Inner Galaxy's GeV excess. Each of these arguments strongly disfavor 
millisecond pulsars as the source of this signal. 

PACS numbers: 97.60.Gb, 95.55.Ka, 95.35.+d; FERMILAB-PUB-13-129-A 



I. INTRODUCTION 

Over the past few years, a number of groups using data 
from the Fermi Gamma-Ray Space Telescope have iden- 
tified a spatially extended component of gamma rays in 
the region of the Galactic Center, peaking at energies 
of ^2-3 GeV [JJ-[5]. This signal has been interpreted as 
possible evidence of annihilating dark matter particles, 
and can be well fit by 7-12 GeV particles annihilating to 
t+t~ (possibly among other leptons) [TH3] or by 22-45 
GeV particles annihilating to quarks [IH31 [5] (alterna- 
tively, see Ref. [5]). In either case, the morphology of 
the gamma-ray signal requires a dark matter distribu- 
tion which scales approximately as p oc r~ 1,2 to r~ 1A , 
where r is the distance to the Galactic Center [TJ [5] , and 
an annihilation cross section that is on the order of a 
few times av ~ 2 x 1CU 27 cm 3 /s [1] (up to overall un- 
certainties in the normalization of the halo profile). This 
required dark matter distribution is in good agreement 
with expectations based on hydrodynamical simulations 
(see Ref. [7] and references therein) and is consistent with 
current observational constraints 8 . Similarly, the re- 
quired annihilation cross section is comfortably within 
the range expected for a thermal relic of the Big Bang. 

Until recently, most opinions regarding the Galactic 
Center's GeV gamma-ray signal fell within one of three 
broad categories. First were those who found the dark 
matter interpretation to be compelling, and considered 
the gamma-ray signal to be difficult to explain with as- 



trophysical sources or backgrounds. Second were those 
who argued that a population of a few thousand millisec- 
ond pulsars concentrated very densely within the inner 
tens of parsecs of the Galaxy were at least as likely to 
be responsible for the observed gamma-ray emission (as 
discussed in Refs. [TH3J [S] and advocated for in Ref. [TO]). 
And third, is the not uncommon view that the region of 
the Galactic Center is too astrophysically complex to reli- 
ably isolate or identify any gamma-ray signal that might 
result from annihilating dark matter particles. 

Very recently, however, this situation has changed sub- 
stantially. In Ref. [11J . an analysis of Fermi data was 
presented which identified a component of gamma rays 
from the inner kiloparsecs of the Milky Way, exhibiting a 
spectrum and morphology consistent with the signal pre- 
viously observed from the Galactic Center. But whereas 
earlier studies were able to confidently identify this sig- 
nal only within a few hundred parsecs of the Galactic 
Center, the template technique applied in Ref. [TT] was 
able to clearly detect this signal out to at least ^2-3 kpc 
to the north and south, and with evidence of extension 
out to ^5 kpc [11] , This result is important in that it 
does not rely on observations of the complex region of 
the Galactic Center, and is not highly sensitive to uncer- 
tainties in Fermi's low energy point spread function, thus 
circumventing the concerns of the third group mentioned 
in the previous paragraph. 

In this paper, we discuss whether this GeV excess ob- 
served from the Inner Galaxy [TJJ could potentially be 
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FIG. 1: The period (P) and its rate of change (dP/dt), for pulsars described in the ATNF catalog |12j . Open circles denote 
binary pulsars (most of which are millisecond pulsars), while squares and X's represent radio quiet pulsars and anomalous 
X-ray and/or soft gamma-ray emitters, respectively. All other pulsars are shown as dots. Also shown are contours of constant 
spin-down power (red dashes), magnetic field (blue dots), and characteristic age, r = P/2P (solid black), calculated assuming 
a neutron star mass of 1.4M and a radius of 10 km. See text for discussion. 



explained by a population of gamma-ray pulsars. Here, 
and throughout this paper, we use the phrase "Inner 
Galaxy" to refer to the region within several kiloparsecs 
of the Galactic Center, as opposed to the much smaller 
region (^100 pc) referred to as the "Galactic Center". 
After briefly reviewing some of the most relevant aspects 
of pulsars in Sec. [TTJ we turn our attention in Sec. |III| 
to the gamma-ray spectra observed from known millisec- 
ond pulsars. We find that the spectrum observed from 
these sources is not compatible with that of the Inner 
Galaxy's GeV excess. Independently of spectral argu- 
ments, we demonstrate in Sec. |IV| that the Inner Galaxy's 
GeV excess cannot be accounted for with millisecond 
pulsars without significantly overpredicting the number 
of pulsars that Fermi would have resolved as individual 
point sources. Pulsar distributions which are consistent 
with Fermi's source catalog can account for no more than 
^10% of the observed GeV excess. From either of these 
arguments, we conclude that gamma-ray pulsars cannot 
account for a significant fraction of the Inner Galaxy's 
GeV excess. In Sec. |VJ we briefly summarize our results 
and their implications. 

II. A BRIEF REVIEW OF ORDINARY AND 
MILLISECOND PULSARS 

Pulsars are rapidly spinning neutron stars which 
steadily convert their rotational kinetic energy into ra- 
diation, including potentially observable emission at ra- 



dio and gamma-ray wavelengths. When initially formed, 
pulsars typically exhibit rotational periods on the or- 
der of tens or hundreds of milliseconds, and mag- 
netic field strengths of -10 n -10 13 G (see Fig. [if. 
As a result of magnetic-dipole braking, a pulsar's 
period will slow down at a rate given by P = 
3.3 x 10- 15 (B/10 12 G) 2 (P/0.3s)- 1 , corresponding to 
an energy loss rate of E = 4tt 2 IP/P 3 — 4.8 x 
10 33 erg/s(P/10 12 G) 2 (P/0.3s)- 4 (//10 45 gcm 2 ). As a 
result of this rotational slowing, pulsars steadily become 
less luminous. For very young pulsars, this occurs very 
rapidly. Within a few hundred thousand years, the Crab 
and Vela pulsars will become a thousand times fainter 
than they are at present. After ^10-100 million years, 
such objects slow to a point at which they are no longer 
able to generate radio emission (crossing what is known 
as the pulsar "death line" ) . 

To those pulsars that are gravitationally bound to a 
binary companion, another stage of evolution is possi- 
ble. If at some point in time (likely well after the pulsar 
has lost most of its rotational kinetic energy and become 
faint) the companion enters a red giant phase, accretion 
onto the pulsar and the corresponding transfer of angu- 
lar momentum can dramatically increase the rotational 
speed of the pulsar (to P^ 1.5-100 ms), while also dra- 
matically reducing the magnetic field (to P~10 8 -10 9 G). 
Such millisecond pulsars (MSPs) (also known as spun-up 
or recycled pulsars) can be as luminous as ordinary pul- 
sars, but evolve much more slowly, remaining bright for 
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billions of years. 

Another important difference between ordinary pulsars 
and MSPs is found in their velocity distributions. In 
the process of their formation, pulsars receive substantial 
kick velocities, resulting from small asymmetries in their 
collapse. The average velocity observed among young 
pulsars is ^400 ± 40 km/s [15] . which is much higher 
than is observed among other stellar populations. MSPs, 
however, necessarily consist of neutron stars that have ei- 
ther retained or captured a binary companion, and thus 
must have had unusually weak kick velocities. Further- 
more, the additional mass of the companion star further 
reduces the velocities acquired by these systems. In Ap- 
pendix]^] we discuss this in more detail and estimate that 
instead of average kick velocities of ~400 km/s, MSPs 
should receive average velocities of 10-50 km/s, with the 
precise value depending on the details of the stellar dis- 
tribution being considered. This conclusion is supported 
by the much lower velocities that are observed among 
MSPs [13] , and by the fact that most observed MSPs re- 
side within globular clusters, which have escape velocities 
on the order of only tens of km/s. 

MSPs are of particular interest in this study for two 
reasons. First, the morphology of the gamma-ray signal 
from the Galactic Center is highly concentrated, much 
more so than is expected from the overall stellar dis- 
tribution. In particular, if any stellar population is to 
produce this signal, its members must preferentially be 
located very centrally around the inner tens of parsecs 
surrounding the Galactic Center, with a number density 
that scales approximately as n(r) oc r~ 2 - 4 to r -2 ' 8 [HE]- 1 
This extremely steep distribution would be very difficult 
to accommodate with ordinary pulsars, whose kick veloc- 
ities are more than sufficient to expel the overwhelming 
majority of pulsars from the Galactic Center. With their 
much weaker kick velocities, however, MSPs which form 
in the Galactic Center could potentially remain concen- 
trated in this region. Furthermore, the large numbers 
of MSPs present in globular clusters suggest that they 
are produced in part as a result of stellar encounters (see 
Sec. 3.3 of Ref. [13], and reference therein). If this conclu- 
sion also applies to the Galactic Center, it could explain 
why the morphology of the observed gamma-ray signal 
is so much more centrally concentrated than the over- 
all stellar distribution. Second, the gamma-ray emission 
identified in Ref. [TT] extends to at least ~3 kpc north and 
south of the Galactic Center. By the time that a newly 
formed ordinary pulsar could travel more than a few hun- 
dred parsecs from the location of its birth (likely near the 
Galactic Plane), it will have lost the vast majority of its 
initial rotational energy and become too faint to signif- 



Note that for dark matter, the annihilation rate per volume is 
proportional to the square of the dark matter density, and thus 
the observed morphology of the Galactic Center's gamma-ray 
signal favors pdm <x r -13 , rather than the n oc r~ 2,6 required 
of a stellar distribution. 
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FIG. 2: The combined gamma-ray spectrum from 37 mil- 
lisecond pulsars observed by Fermi. The solid line shows 
the best-fit parametrization to this spectrum, dN-^/dE^ oc 
S7 146 exp(-£ 7 /3.3GeV). 



icantly contribute to the observed gamma-ray emission. 
Observations also support the fact that ordinary pulsars 
are found overwhelmingly within or near the volume of 
the stellar disk. In contrast, MSPs can remain bright for 
billions of years. Thus if a significant number of MSPs 
somehow escape the gravitational potential of their en- 
vironment, there could potentially exist a population of 
luminous, high-latitude gamma-ray pulsars. 



III. COMPARISON WITH THE OBSERVED 
GAMMA-RAY SPECTRA OF MILLISECOND 
PULSARS 

The Fermi Collaboration has detected gamma-ray 
emission from a total of 125 sources identified as pul- 
sars, 47 of which have millisecond-scale periods [T3]. Of 
these 47 MSPs, 37 have spectral information listed in the 
second Fermi source catalog (2FGL) [T7]. In Fig. [2j we 
plot the sum of the spectra of these 37 sources. We find 
that this collection of sources can be very well fit by the 
standard pulsar spectral parametrization, dN^/dE^ oc 
E~ T exp(— E 1 /E cut ), with best-fit values of T = 1.46 and 
E C ut — 3.3 GeV. This is very similar to the results found 
in earlier studies, based on a smaller number of Fermi 

msps pa Eg. 

In Fig. [3| we compare this best-fit spectrum (solid line) 
to that observed from the |6| = 10° — 20° regions of 
the Inner Galaxy, after subtracting an inverse Compton 
component that accounts for the Fermi Bubbles emission 
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FIG. 3: The gamma-ray spectrum of the |6| = 10° — 20° re- 
gions of the Inner Galaxy (see Ref. [II]), after subtracting 
emission from inverse Compton scattering compared to 
the spectral shape best-fit to 37 MSPs observed by Fermi 
(solid line); see Fig. [5] Also shown for comparison are 
the shapes corresponding to spectral parameters that better 
match this emission: V — 1.0, E cut = 2.75 GeV (dashed) and 
T = 0.5, Scut = 2.0 GeV (dotted). 



(see Ref. [UJ for details). 2 The spectral shape observed 
from these 37 resolved MSPs exhibits a much softer spec- 
tral index than the spectrum of the excess emission ob- 
served from the Inner Galaxy, especially at energies be- 
low ~l-2 GeV. Also shown for comparison are harder 
spectral shapes, corresponding to T = 1.0, E cut = 2.75 
GeV (dashes) and T = 0.5, E cut = 2.0 GeV (dots). 
While such hard spectra can provide a good fit to the 
emission observed from the Inner Galaxy (especially the 
T = 0.5, E cut = 2.0 GeV case), they are not consistent 
with the spectral shape shown in Fig. [2j The compar- 
ison between these harder spectral shapes and the er- 
ror bars shown in Fig. [2] yields fits of % 2 = 17.8 and 
38.9 (over 5-1 degrees-of-freedom) for these two parame- 
ter sets (T = 0.5, E cut = 2 GeV, and T = 1.0, E cut = 2.75 
GeV, respectively), each of which can be excluded at be- 
yond the 99.8% confidence level. At least at face value, 
it appears that we can exclude at high confidence a MSP 
origin for the emission observed from the low-latitude re- 
gions of the Inner Galaxy. 

Perhaps, however, the MSPs that have been resolved 
by Fermi are not representative of all such objects, and 
the diffuse emission from the sum of all unresolved (and 



2 Note that while we have chosen to compare to the spectrum 
observed from the |fe| = 10° — 20° region of the Inner Galaxy, this 
spectrum is very similar to that extracted from higher latitude 
regions |11| . We have chosen to not use the spectrum extracted 
from the |b| = 1° — 10° region due to difficulties in separating 
this signal from emission associated with the Galactic Disk. 



faint) MSPs has a much harder spectral index than is 
observed from resolved (and brighter) sources. Among 
Fermi's MSPs, however, we see no evidence for this. In 
particular, we find a best-fit spectral index of T = 1.36 
for the sum of the 21 MSPs with gamma-ray luminosi- 
ties less than 10 37 GeV/s (above 100 MeV), and V = 1.34 
for the sum of the 8 MSPs with gamma-ray luminosities 
less than 10 36 GeV/s. These values are not very differ- 
ent from that found in our overall best- fit, T = 1.46. If 
MSPs exhibit a correlation between hard spectral indices 
and low luminosities, this trend is not evident among the 
observed source population. 

Furthermore, if hard spectrum, low-luminosity sources 
dominate the diffuse emission from MSPs, then the hard 
spectral index should be reflected in the emission ob- 
served from globular clusters, which should contain a 
representative sample of MSPs. Although the spectra of 
the 11 globular clusters included in the 2FGL [17] have 
large error bars and are thus difficult to evaluate indi- 
vidually, the sum of the spectra from these 11 sources is 
quite similar to that observed from Fermi's MSPs. Simi- 
larly, the Fermi Collaboration studied 8 globular clusters 
and found their (statistically weighted) average spectral 
index to be T =1.35 [T5], again similar to that observed 
from resolved MSPs. 

Although we have shown in this section that the 
gamma-ray spectrum observed from individual MSPs 
(and from globular clusters) is incompatible with that 
from the Inner Galaxy as reported in Ref. [11] , one might 
worry that systematic uncertainties in the low-energy 
(< 1 GeV) spectrum could possibly alter this conclusion. 
The error bars presented in Ref. |llj (and shown in our 
Fig. [3]) are purely statistical, and do not reflect the pos- 
sible mismodeling of point sources or of diffuse emission. 
While the over-subtraction of low-energy emission from 
known point sources could, in principle, lead to an ar- 
tifically hard spectrum at low-energies, only if the Fermi 
collaboration's source catalog [2U] overestimates the total 
flux from the 35 sources in the \b\ — 10° — 20° region, for 
example, by more than a factor of two in the in the 300- 
1000 MeV range (a variation several times larger than the 
quoted errors) could the spectrum of the Inner Galaxy's 
GeV excess be consistent with that observed from indi- 
vidual MSPs. More difficult to rule out is the possibily 
that the Fermi collaboration's diffuse model significantly 
overestimates the density of cosmic rays in the region of 
interest, leading the analysis of Ref. [TT] to effectively 
oversubtract gamma-ray emission from pion production 
and other diffuse processes, potentially artificially hard- 
ening the spectrum of the GeV excess at low energies. We 
note that the fit residuals from the analysis of Ref. [TT] , 
averaged over the regions in question, are much smaller 
than the signal at all relevant energies; re-adding them 
to the signal does not meaningfully soften the spectrum. 

Although systematic uncertainties in the Fermi instru- 
ment response functions below 1 GeV could plausibly 
skew the inferred spectral shape, no evidence for this is 
seen in other spectral components, such as that associ- 
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FIG. 4: The observed flux distribution (proportional to dN/d\ogS) of identified millisecond pulsars with |6| > 10° (solid black), 
compared to that predicted in the base model of Ref. [2T] ((|z|) = 1 kpc, ay = 5 kpc, Bo = 10 8 G, normalized to accomodate 
the observed number of very bright sources). Also shown are the distributions of identified MSPs plus all unidentified Fermi 
sources (dotted blue), and of identified MSPs plus all unidentified sources found by the Sibyl algorithm [27] to be either likely 
pulsars or sources of an inconclusive nature (dashed black). Also shown is the range of Fermi's threshold to resolve an individual 
source [17] . This base model cannot account for the observed number of bright MSPs without significantly overpredicting the 
number of fainter MSPs. 



ated with the Fermi diffuse model [TTJ ; this argues against 
an energy-dependent error in Fermi's effective area be- 
ing responsible for the apparently hard spectrum of the 
Inner Galaxy's GeV excess. Furthermore, the large an- 
gular size of the regions of interest, and their significant 
distance from the Galactic Center, make any mismodel- 
ing of Fermi's point spread function an unlikely source of 
large distortions to the spectrum. 

To summarize the results of this section, we find that 
the gamma-ray spectra observed from individual MSPs 
consistently reveal a spectral index that is much too 
soft to accommodate the signal observed from the Inner 
Galaxy. Furthermore, we find no evidence for a popu- 
lation of low-luminosity and spectrally hard MSPs that 
might be able to account for the signal. 



IV. THE DISTRIBUTION OF MILLISECOND 
PULSARS IN THE MILKY WAY 

In the previous section, we showed that the gamma- 
ray spectrum observed from individual MSPs (and from 
collections of MSPs in globular clusters) is not consistent 
with the spectral shape of the Inner Galaxy's GeV excess. 
In this section, we set aside this conclusion for the time 
being and focus instead on constraints derived from the 
observed spatial and flux distributions of MSPs. We will 



use this information to assess the question of whether 
the intensity and morphology of the Inner Galaxy's GeV 
excess might originate from a population of unresolved 
MSPs. 



A. Millisecond Pulsars Associated with the 
Galactic Disk 

We begin by considering MSPs which follow a distri- 
bution similar to that of the Milky Way's disk. As our 
starting point, we adopt the "base model" of Ref. |21) . 
which includes a spatial distribution and luminosity func- 
tion for MSPs in the Milky Way. In particular, we adopt 
a spatial distribution of MSPs with a number density 
given by: 

n(r,z) cx exp(-r 2 /2^) exp(-\z\/{\z\)), (1) 

where r and z describe the location in cylindrical coor- 
dinates. To begin, we will consider values of a r = 5 kpc 
and (|z|) = 1 kpc, as adopted in the "base model" of 
Ref. [21] • 

Again following Ref. [21], we take the gamma-ray lu- 
minosity (above 100 MeV) of a MSP to be equal to 5% 
of its energy loss rate, E, except for the most luminous 
sources which follow cx \pE~. For the distribution of 
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MSP periods, we assume dN/dP oc P~ 2 , with a mini- 
mum value of 1.5 msec [22] (the most rapidly spinning 
pulsar observed to date has a period of 1.4 msec [23] ) - 
The time derivative of a MSP's period is determined by 
its magnetic field (through magnetic dipole braking, see 
Sec. The magnetic fields are taken to follow a log- 
normal distribution centered around B = 10 s G and 
with a logarithmic standard deviation of 0.2. While we 
take the gamma-ray spectrum to follow the form of the 
best-fit as shown in Fig. [2j the precise spectral shape 
does not significantly impact any of the results presented 
in this section. 

In Fig. |4j we show the flux distribution of high-latitude 
(|6| > 10°) MSPs (proportional to dN/d\ogS) pre- 
dicted in the base model of Ref. 21 (labeled "FGL Base 
Model"). We compare this prediction to the number of 
sources as observed by Fermi. Here, the solid black his- 
togram describes the distribution of sources in the 2FGL 
which have been identified as MSPs [16] , while the dotted 
blue histogram denotes the sum of the identified MSPs 
along with all presently unidentified sources in the cat- 
alog (i.e. all sources listed as unassociated in the 2FGL 
that do not appear on the list of Fermi pulsars [16] and 
have not since been identified in Ref. [23] as a blazar) . For 
the predicted distribution, we have normalized the total 
number of MSPs to approximately match the observed 
number of very bright MSPs (F 7 (> 1 GeV) ~ 1(T 8 ph 
cm~ 2 s _1 ). Note that in this respect, we depart from the 
base model of Ref. [21] . For this choice of normaliza- 
tion, we find that unresolved MSPs in this model pro- 
duce about 0.5% of the high latitude (|6| > 40°) diffuse 
gamma-ray background at E~ ~ 1 GeV, which is about 
a factor of three below the maximum value consistent 
with Fermi's anisotropy constraint [35] [2S]. The verti- 
cal dashed lines in Fig. [4] denote the range of Fermi's 
threshold for a source out of the plane (|6| > 10°) to 
be included in the 2FGL catalog, as quoted in Ref. [T7] , 
The range of this threshold spans sources with effective 
spectral indices between -1 and -3 (assuming a power-law 
form) . We note that the point source threshold for MSPs 
should typicaly fall near the lower end of this range, since 
their spectra are relatively hard at ~1 GeV. 

While this model predicts a diffuse gamma-ray signal 
from the Inner Galaxy that is similar to the observed GeV 
excess (with a similar morphology, and an overall inten- 
sity that is only a factor of a few less than observed), it 
also significantly overpredicts the number of MSPs with 
F 7 (> lGeV) - 10~ 9 ph cm^s" 1 . Only if essentially 
all of Fermi's unidentified sources (above |6| = 10°) are 
MSPs could this model be potentially compatible with 
the observed flux distribution. It is clear, however, that 
only a modest fraction of these unidentified sources are 
pulsars, and that most of them are instead blazars or 
other types of active galactic nuclei (AGN). For exam- 
ple, the authors of Ref. [27], using the random forest 
classifier Sibyl, trained on the observed spectra and vari- 
ability of over 900 identified Fermi point sources (AGN 
and pulsars), determined that at least 80% of Fermi's 
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FIG. 5: The combined gamma-ray spectrum from all of 
Fermi's unidentified sources located outside of the plane (|6| > 
10°). The solid line shows the best-fit parametrization to this 
spectrum, dN 1 /dE 1 oc E~ 2 ' 22 exp(— _E 7 /33 GeV). This spec- 
trum does not resemble that observed from individual mil- 
lisecond pulsars, but is consistent with that observed from 
blazars and other types of AGN. 



unidentified high-latitude > 10°) sources are likely 
AGN. Furthermore, the overall spectrum from this col- 
lection of unidentified sources does not resemble that ob- 
served from individual MSPs (or observed from the Inner 
Galaxy), but instead resembles that of AGN. In Fig.[5j we 
plot the combined spectrum of these unidentified sources 
(all with |6| > 10°). This spectrum shows no sign of a 
sharp spectral peak at ~2 GeV, as the Inner Galaxy's dif- 
fuse emission does, nor does it resemble the more mildly 
peaked spectrum observed among the identified MSPs. 
The shape of this spectrum strongly suggests that most 
of the unidentified sources are not pulsars, but are instead 
mostly AGN or other soft-spectrum gamma-ray sources. 

In Fig. [4] the black-dashed histogram represents the 
distribution of identified MSPs added to the distribution 
of sources classified by Sibyl as either a likely pulsar, or 
as a source of an inconclusive nature (only sources classi- 
fied as likely AGN were not included in this distribution) . 
This distribution represents an approximate upper limit 
for the numbers of Fermi's sources that could potentially 
be MSPs. In all likelihood, the true distribution falls 
somewhere between the solid-black histogram (presently 
identified MSPs) and the dashed-black histogram (iden- 
tified MSPs plus Sibyl's likely pulsars and inconclusive 
sources). When comparing these distributions to that 
predicted by the base model of Ref. [2T] , we are forced to 
conclude that this model cannot account for the observed 
number of very bright MSPs without predicting far too 
many fainter MSPs. 

To better accommodate the observed flux distribution 
of MSPs, we must consider population models with ei- 
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FIG. 6: As in Fig. [4] but for a number of variations in the millisecond pulsar population model. See text for details. 



ther 1) spatial distributions which are more weighted 
more toward nearby MSPs, or 2) luminosity functions 
that are more weighted more toward higher luminosity 
MSPs. In Fig. |6j we show how varying a number of our 
model's parameters can impact the flux distribution of 
MSPs. In the upper frames, we vary the parameters of 
our spatial distribution, (\z\) and oy (see Eq. [I]). By re- 
ducing the vertical scale height of the MSPs distribution 
to (\z\) = 0.3 kpc (approximately the scale height of the 
Milky Way's thin disk) , the model can provide a not un- 
reasonable match to the observed distribution (although 
nearly all of Sybil's non-AGN sources would have to be 
MSPs in this case). Values of (\z\) > 1 kpc appear to 
be incompatible with the observed flux distribution. In 
contrast, reasonable variations in a r have relatively lit- 
tle impact on the predicted distribution. In the lower 
left frame, we consider the possibility of a local overden- 
sity of MSPs (enhanced by a factor of 10 within 0.3 kpc 
of the Solar System). This, however, had little impact 
on the overall distribution, except for slightly increasing 



the predicted number of very bright sources. Lastly, in 
the lower right frame we focus on the MSP luminosity 
function by varying in the central value of the magnetic 
field distribution, Bq. We find that by increasing this 
quantity from 10 8 to 10 8 ' 5 gauss or higher, we can much 
better accommodate the observed flux distribution. We 
also note that from the information shown in Fig. [I] val- 
ues of B ~ 10 8 3 — 10 8 5 G appear to best describe the 
observed population of MSPs. 

The spatial distribution of MSPs is not entirely un- 
constrained, however. In Figs. [7] and [5J we plot some of 
the information we have pertaining to the spatial distri- 
bution of the MSPs observed at radio and gamma-ray 
wavelengths, respectively. Fig. [7] shows the distance to 
pulsars, and the distance of those pulsars from the Galac- 
tic Plane, as a function of period, for pulsars in the ATNF 
catalog [12]. The collections of points forming horizon- 
tal lines in these plots are groups of millisecond pulsars 
found in globular clusters. In Fig. [8] the spatial distri- 
bution of those MSPs observed by Fermi and with coor- 
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FIG. 7: The distance from the Solar System (left) and distance from the Galactic Plane (right), as a function of period, for 
pulsars in the ATNF catalog [12] . The groups of points that form horizontal lines are millisecond pulsars in globular clusters. 
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FIG. 8: Left: The distance from the Galactic Plane and distance to the Solar System projected along the Galactic Plane for 
those millisecond pulsars observed by Fermi and with coordinates given in the ATNF catalog. Right: The locations in the 
Galactic Plane of those millisecond pulsars observed by Fermi and with coordinates given in the ATNF catalog, and with 
\z\ < 1.5 kpc. In this coordinate system, the Galactic Center is located at (0,0), while the Solar System is at (0,8.5 kpc). 
Shown for comparison are the approximate locations of the Orion-Cygnus, Sagittarius, and Perseus arms of the Milky Way. 



dinates given in the ATNF catalog is shown. From these 
two figures, it is clear that the MSP distribution is not 
highly concentrated within a few hundred parsecs of the 
Galactic Plane, instead favoring a distribution at least 
as broad as (|z|) ~ 0.5 kpc. 3 Thus reducing the value of 
alone (as shown in the upper left frame of Fig. [6]) 
does not seem to be a viable way to accommodate the 
observed flux distribution. Furthermore, the left frame 



3 While (|z|) ~ 0.5 kpc provides the best-fit to the observed dis- 
tribution of observed MSPs, observational bias favoring nearby 
sources might lead us to slightly underestimate this quantity. We 
take 0.5 kpc to be an approximate lower limit for 



of Fig. [7] and the lower frame of Fig. [8] do not reveal any 
very large local overdensity of MSPs. 

From the results shown in Fig. [6j constrained by the 
observed spatial distributions of Figs. [jj and [HJ we con- 
clude that we are forced to consider MSP luminosity func- 
tions favoring somewhat higher values than are found in 
the FGL base model. In terms of the magnetic field pa- 
rameter, this favors Bq ~ (2 — 6) x 10 s G, although one 
should keep in mind that this parameter is somewhat de- 
generate with the period distribution, and with the frac- 
tion of rotational energy loss that goes into gamma ray 
production. We consider examples of what appear to be 
reasonably viable MSP population models in Fig. [9j In 
the upper frames, we show the flux distribution; for the 
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FIG. 9: Top: As in Fig. [6j but for parameters which yield flux distributions which are in reasonable agreement with observations. 
Bottom: The observed gamma-ray flux (after subtracting inverse Compton emission) between 1.9 and 3.5 GeV from the regions 
associated with the Fermi Bubbles, in five latitude bands (jbj = 1° - 10°, 10° - 20°, 20° - 30°, 30° - 40°, and 40° - 50°), 
compared to the prediction from MSPs in the same four models used in the upper frames. In each case, only ~5-10% of the 
observed emission can be accounted for by millisecond pulsars. See text for details. 



lower choice of Bq used in each frame (dot-dashed), we 
approximately saturate the observed source distribution. 
The distributions shown for slightly larger values of Bq 
should be considered more realistic, many of the sources 
included in the dashed histogram are likely to be sources 
other than MSPs (in particular among Sybil's inconclu- 
sive sources). In the lower frames of Fig. [9j we show the 
gamma-ray flux at 1.9-3.2 GeV (the approximate peak of 
the observed excess) observed by Fermi from various lat- 
itude ranges of the Inner Galaxy [llj . and compare this 
to the predicted flux in the same four MSP population 
models. Clearly these models cannot account for the ob- 
served emission, falling short in each case by a factor of 
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4 In calculating the contribution to the diffuse gamma-ray emis- 
sion as shown in the lower frames of Fig. [9] we have treated 
any MSP with a flux less than 4.1 X 10~ 10 cm" 2 s" 1 above 1 
GeV as unresolved and included its emission in the prediction 
for the diffuse flux. In light of the hard spectra of MSPs, this is 
a fairly conservative threshold, and we expect most MSPs with 
fluxes above ~ 2.4 X 10 — 10 cm -2 s — 1 to be resolved (see Fig. 6 
of Ref. 17 ). At high-latitudes, this provides a reasonable upper 
limit for the contribution to the diffuse flux. At lower-latitudes, 
however, some MSPs slightly brighter that our assumed thresh- 
old may go unresolved. If we increase our point source threshold 
by a factor of 2 (as is appropriate for sources at |6| ~ 10° |20|). 
we find that the low-latitude diffuse flux approximately doubles, 
still falling well short of that required to explain the observed 



10 




FIG. 10: As in Figs. [4] [6] and [9] but now also showing the contribution from a population of millisecond pulsars associated with 
the Galactic Bulge. The three solid blue lines correspond to spatial distributions which are a spherical gaussian with ct_r = 0.5, 
1 and 3 kpc (from bottom-to-top, although the cr^j = 0.5 contour falls below the range shown in the right frame). We have 
normalized the bulge contribution such that the number of millisecond pulsars per stellar mass is the same in the bulge as in 
disk. Here, we have also adopted a disk distribution with (\z\) = 0.5 kpc. 



B. Millisecond Pulsars Associated with the 
Galactic Bulge 

In the previous subsection, we showed that after tak- 
ing into account the observed flux distribution and spa- 
tial distribution of MSPs, the population of MSPs asso- 
ciated with the Galactic Disk cannot produce more than 
~5-10% of the Inner Galaxy's GeV excess. To increase 
the intensity of diffuse emission from unresolved MSPs 
in the Inner Galaxy without predicting far too many 
resolved sources, we must require an additional popu- 
lation of MSPs with a flux distribution that increases 
sharply below the point source sensitivity of Fermi. The 
most promising way to accomplish this is to add an ad- 
ditional component to our population model which ex- 
plicitly takes into account those MSPs associated with 
the Galactic Bulge. In this subsection, we consider MSP 
models which include contributions from sources associ- 
ated with both the disk and the bulge of the Milky Way. 

We model the bulge population of MSPs with the same 
luminosity function as the disk component, and with 
a spatial distribution that is described by a spherically 
symmetric gaussian, n(R) oc exp(—R 2 /a Fl ), where R is 
the distance to the Galactic Center (or is not to be con- 
fused with the quantify ay, as appears in Eq. [lj. In 
Fig. [lO] we show the flux distribution for disk and bulge 
components, using three values of ctr, and for B = 10 8 ' 25 
G (left), B = 10 8 ' 5 G (right), and (|z|) = 0.5 kpc. If we 
adopt a MSP distribution for the bulge that is similar 



emission. 



to the distribution of bulge stars (<tr ~ 0.5 kpc), we get 
almost no contribution (the bottom blue curve barely ap- 
pears in the left frame and falls below the range shown 
in the right frame, and thus does not appear) . If we in- 
crease ctr to 1 kpc or more, we find a significantly larger- 
contribution from the bulge, but also a non-negligible 
contribution to the number of individual sources that 
should be resolvable by Fermi. In particular, three of 
Fermi's observed MSPs exhibit gamma-ray luminosities 
of 2 x 10 37 GeV/s or higher, each of which would be well 
above Fermi's point source threshold if located at a dis- 
tance of ~10 kpc from the Solar System (for \b\ > 10°). 
The fact that a non-negligible fraction of bulge MSPs 
(those with high luminosities and outside of of the Galac- 
tic Plane) will be resolvable as individual point sources 
by Fermi will ultimately limit how much of the Inner 
Galaxy's GeV excess we can attribute to such a popula- 
tion. 

The bulge contributions shown in Fig. [10] have been 
normalized assuming that the number of MSPs per stel- 
lar mass is the same in the bulge as in the disk. It is 
possible, however, that the relative number of MSPs in 
the bulge could be somewhat larger. As an extreme illus- 
tration of this possibility, we note that globular clusters 
are observed to contain ~10 2 times more low mass X-ray 
binaries (the assumed progenitors of MSPs) per stellar 
mass than is found throughout the disk [28, 29 . This 
is presumably the consequence of the very high stellar 
densities found in these systems (up to ~10 2 -10 3 stars 
per cubic parsec, compared to ^0.4 in the local volume 
of the disk), which can be expected to significantly in- 
crease the probability that a given pulsar will obtain a 
companion and thus potentially evolve into a MSP. The 
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FIG. 11: Top: As in Fig. |10| but now also showing the bulge, disk, and bulge+disk contributions from millisecond pulsars. 
Here, we have adopted an =1 kpc and = 0.5 kpc. We have normalized the bulge contribution such that the number of 
millisecond pulsars per stellar mass is the same in the bulge as in disk (solid blue and solid red) and by a factor that is 2 times 
larger (dashed red). Bottom: As in the lower frames of Fig. [9] but for the sum of disk and bulge contributions. The total 
diffuse emission from millisecond pulsars is in each case found to be much less than that needed to account for the observed 
GeV excess. 



average stellar density in the bulge is significantly higher 
than in the disk, but much lower than that found in the 
cores of globular clusters (only in the innermost tens of 
parsecs around the Galactic Center is the stellar density 
comparable to that found in globular clusters). As a re- 
sult, we naively expect only a modest enhancement in 
the number of MSP per stellar mass found in the bulge 
relative to that in the disk (likely on the order of a few 
or less) [9l [30] . This conclusion is further supported by 
the observed distribution of low mass X-ray binaries in 
the Galactic Bulge [5T] , 



In Fig. 11 we show the distribution of sources (top) and 
flux of diffuse gamma-ray emission (as a function of lati- 
tude) from MSPs, including contributions from the disk 
and bulge. Here we have chosen a bulge distribution de- 



scribed by Or = 1 kpc because significantly smaller val- 
ues lead to a negligible contribution to the diffuse emis- 
sion, while much larger values predict numbers of ~10 
cm~ 2 s _1 sources that exceed those observed by Fermi. 
Furthermore, we find that our conclusions are not sensi- 
tive to the precise value of this parameter. We have also 
taken here a disk width of (|z|) = 0.5 kpc, which approxi- 
mately maximizes the allowed contribution to the diffuse 
emission from MSPs in the bulge. In the lower frames of 
this figure, we show the gamma-ray flux at 1.9-3.2 GeV 
(the approximate peak of the observed excess) observed 
by Fermi from various latitude ranges of the Fermi bub- 
bles |llj , and compare this to the predicted flux in these 
disk+bulge MSP population models. Again, we find that 
the predicted contribution from MSPs cannot account for 
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the observed emission, falling short in each case by a fac- 
tor of -10. 



V. SUMMARY AND CONCLUSIONS 

Millisecond pulsars have been proposed as a possi- 
ble source for the Inner Galaxy's GeV excess identified 
within the data of the Fermi Gamma-Ray Space Tele- 
scope. This hypothesis has been motivated by two main 
considerations. First, the spectrum of gamma-ray pul- 
sars is observed to peak at E 1 ~ 1-2 GeV, similar to 
that of the GeV excess. Second, it is plausible that the 
high stellar densities found in the Galactic Center could 
facilitate the production of a large number (~ 10 3 ) of mil- 
lisecond pulsars, with a spatial distribution that is very 
highly concentrated within the innermost tens of parsecs 
of the Milky Way. It is far less clear, however, that such 
objects could account for the angular extent of this ex- 
cess, which has recently been shown to extend out to at 
least ~3 kpc from the Galactic Center [TT]. In this pa- 
per, we address specifically this more extended signal, 
and the question of whether pulsars might account for 
this observed emission. 

In considering the possibility that a large population of 
unresolved millisecond pulsars is the source of the Inner 
Galaxy's GeV excess, we have presented two independent 
arguments, each of which strongly disfavors this hypoth- 
esis: 



• In Sec. Ill we showed that the spectrum of the 
millisecond pulsars resolved by Fermi is not com- 
patible with the observed spectral shape of the 
GeV excess. In particular, the combined spec- 
trum of Fermi's 37 millisecond pulsars with spec- 
tral information contained in the 2FGL catalog is 
well fit by dN^/dE 1 cx E~ 1A6 exp(— i? 7 /3.3 GeV), 
while the spectrum of the GeV excess is much 
harder at sub-GeV energies, i.e. dN 7 /dE 7 cx 
E~ - 5 exp(- J B 7 /2.75 GeV). 



• In Sec. |IV[ we considered models to describe the 
spatial distribution and luminosity function of mil- 
lisecond pulsars in the Milky Way. After con- 
sidering a wide range of parameters and distribu- 



tions, we found no models that could accommodate 
the GeV excess without significantly overpredict- 
ing the observed number of bright, high-latitude 
(\b\ > 10°) millisecond pulsars. Models that did 
not violate this constraint were capable of produc- 
ing no more than ~10% of the GeV excess. 

In light of these results, we are forced to conclude that 
millisecond pulsars are not responsible for the GeV ex- 
cess observed from the Inner Galaxy. Although one could 
imagine another (unknown) class of gamma-ray sources 
that could account for this signal, its members would 
be required to have a number of rather specific charac- 
teristics: 1) A very hard spectrum (much harder than 
pulsars), peaking at 2-3 GeV, 2) Low gamma-ray lu- 
minosities (to avoid being identified as individual point 
sources), consistently less than ~10 37 GeV/s, and 3) A 
spatial distribution that is concentrated around the Inner 
Galaxy, with significantly more sources associated with 
the bulge than the disk (per stellar mass). No known 
class of gamma-ray sources exhibits these characteristics 
and, to the best of the authors' knowledge, no such class 
of sources has been proposed. 

In excluding gamma-ray pulsars as the source of the 
Inner Galaxy's GeV excess, the results presented in this 
study provide further support for a dark matter interpre- 
tation of this signal. At this time, we know of no viable 
alternative to dark matter annihilations as the source of 
the excess GeV emission observed from the Galactic Cen- 
ter and from the inner kiloparsecs of the Milky Way. 
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Appendix A: Kick Velocities and Millisecond Pulsars 

Stars massive enough to form a neutron star (M > 
8M©) become supergiants prior to collapse. If such a 
star is to leave its binary companion intact, the stars 
must be separated by a distance of at least an AU 
or so. With this picture in mind, we start by con- 
sidering a binary system with masses Mi and M comp 
and in an orbit with a semi-major axis, a. The initial 
speed of the first star (which will become the pulsar), 
is given by \Vi\ = y/G(Mi + M comp )/a. Upon collaps- 
ing into a neutron star, the first star is given a kick 
velocity, Vkick- Adding these velocities, we arrive at 
V f = V k i C k + Vi. If \V f \ < \Vi\yfi/x, the binary sys- 
tem will remain intact, potentially allowing for the fu- 
ture formation of a millisecond pulsar. The quantity 
X = (Mi + M comp )/(Mf + M comp ) takes into account 
the significant degree of mass loss that is expected to 
occur in explosion (typically Mi ~ 8 — 20 M while 
M f ~ 1.4- 3.2 M Q ). 

As a numerical example, we consider a binary system 
with initial masses of 10 Mq and 1.2 Mq, and in an orbit 
with a semi-major axis of 1 AU. Based on observed pul- 
sar velocities, we take the distribution of initial kicks to 
be isotropic and maxwellian, with a mean velocity of 400 
km/s. For these parameters, and for a final mass of the 
neutron star of Mf = 1.5 M©, we find that only about 
1% of such systems will remain bound to each other (the 
other 99% become ordinary pulsars, unless they later be- 
come bound to a new stellar companion). Of those pul- 
sars which remain bound to their original companion, 
the mean kick velocity of the resulting binary system is 
only 27.6 km/s, which is sufficiently low to explain why 
so many observed MSPs are found within the relatively 
weak gravitational potentials of globular clusters. And 
while the details of this result depend on the input dis- 
tributions of stellar masses, orbits, and initial kick veloc- 
ities, the conclusion that MSPs should typically receive 
kick velocities that are smaller than those of ordinary 
pulsars by an order of magnitude or more is a generic pre- 



diction. In particular, for a = 1 AU and M t 



comp 



1.2 M, 







we find average MSP kicks of 42 km/s for very massive 
neutron stars (Mi = 30 Mq , Mf = 3.2 Mq) and 26 km/s 
for minimally massive (Mi = 8M , Mf = 1.4M ) neu- 
tron stars. Increasing a or M comp further reduces the 
average kick velocity. We do not consider much smaller 
values of M comp than 1.2 M , as such stars would not 
have yet evolved into a giant phase, and thus would not 
have yet spun-up their neutron star companion into a 
MSP. 

We also note that it has been suggested that the distri- 
bution of pulsar kick velocities may be bi-modal. Specif- 
ically, electron-capture supernovae are thought to often 
result in pulsars with relatively weak kick velocities [32] . 
This could significantly enhance the fraction of pulsars 
which retain a binary companion and thus potentially 
evolve into a MSP. 



